Dynamical response of the Hodgkin-Huxley model in the high-input regime 
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The response of the Hodgkin-Huxley neuronal model subjected to stochastic uncorrelated spike 
trains originating from a large number of inhibitory and excitatory post-synaptic potentials is ana- 
lyzed in detail. The model is examined in its three fundamental dynamical regimes: silence, bista- 
bility and repetitive firing. Its response is characterized in terms of statistical indicators (interspike- 
interval distributions and their first moments) as well as of dynamical indicators (autocorrelation 
functions and conditional entropies). In the silent regime, the coexistence of two different coherence 
resonances is revealed: one occurs at quite low noise and is related to the stimulation of subthreshold 
oscillations around the rest state; the second one (at intermediate noise variance) is associated with 
the regularization of the sequence of spikes emitted by the neuron. Bistability in the low noise limit 
can be interpreted in terms of jumping processes across barriers activated by stochastic fluctuations. 
In the repetitive firing regime a maximization of incoherence is observed at finite noise variance. 
Finally, the mechanisms responsible for spike triggering in the various regimes are clearly identified. 

PACS numbers: 05.45.-a,87.10.+e,87.17.Aa,05.40.Ca 
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I. INTRODUCTION 

Neuronal models represent a fundamental benchmark 
to investigate the dynamical response of excitable sys- 
tems under the influence of noise. One of the main rea- 
sons justifying the interest of neuroscientists for this sub- 
ject resides in the observation that in vivo neocortical 
neurons are subjected to a constant bombardment of in- 
hibitory and excitatory post-synaptic potentials (EPSPs 
and IPSPs), somehow resembling a background noise. As 
a consequence in the last decades a large number of nu- 
merical and theoretical studies have been devoted to the 
characterization of the response of simple and more elab- 
orated neuronal models under the influence of a large 
variety of stochastic inputs 0, 0, 13 ■ 

Among the many proposed biophysical models the one 
introduced by Hodgkin and Huxley in 1952 [4| can still 
be considered as a valid framework for exploring neural 
excitability, due to its relative simplicity combined with 
the fact that it embodies the major features of mem- 
brane potential evolution In particular, in order to 
understand the origin of the variability observed in the 
distribution of spikes emitted by cortical neurons the 
response of the Hodgkin-Huxley (HH) model has recently 
been studied under the influence of additive noise B,|p or 
subjected to trains of post-synaptic potentials [7ll9l.ll0|. 

Stochastic resonance (SR) and coherence reso- 
nance (CR) fl2l Il3| represent some of the most inter- 
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esting phenomena observed experimentally and numer- 
ically for excitable neuronal systems driven by noise. 
While SR is related to the enhanced ability of neurons 
to detect weak (periodic or aperiodic) signals when sub- 
jected to additive noise, CR refers to the regularization 
of the response of the system at an optimal noise in- 
tensity in the absence of an external signal (for a com- 
prehensive review see For neuronal systems evi- 
dences of CR have been reported experimentally for the 
cat's spinal and cortical neural ensembles |l5j and the- 
oretically for the following models: FitzHugh-Nagumo 
Fill, leaky integrate-and-fire 0, 0, Hindmarsh-Rose 
18] and Morris-Lecar 111. CR has also been observed 
for the HH model |20l l21j , but these results mainly refer 
to additive continuous noise and to the silent regime near 
the saddle-node bifurcation of limit cycles. 

Our aim is to perform a detailed analysis of the re- 
sponse of the HH model subjected to many stochastic 
trains of EPSPs and IPSPs in its three fundamental dy- 
namical regimes: the silent, the bistable and the repet- 
itive firing ones. The neuron is studied in the so called 
high-input regime @, i.e. when it receives hundreds or 
thousands of post-synaptic inputs per emitted spike. In 
this situation the stochastic input can simply be char- 
acterized in terms of its average, representing the bifur- 
cation parameter of the model, and its variance. Most 
attention is devoted to the mechanisms responsible for 
neuronal firing in the different dynamical regimes and to 
the characterization of CR in terms of dynamical and 
statistical indicators. In particular, we show that the 
conditional entropies can be employed as powerful indi- 
cators to detect coherence, similarly to what has previ- 
ously been done to characterize stochastic resonance for 
a Schmitt trigger . Sound evidences of two coexisting 
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a x {V) (s- 1 ) 


0x(V) (s- 1 ) 
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0.1 (V+40)/(l-exp[-(V+40)/10]) 


4 exp[-(V+65)/18] 


11 


0.01 (V+55)/(l-exp[-(V+55)/10]) 


0.125 exp[-(V+65)/80] 


h 


0.07 exp[-(V+65)/20] 


l/(exp[-(V+35)/10]+l) 



TABLE I: ax{V) and f3x(V) functions (X — m, n, h) entering 
in Eq. @ for the voltage expressed in mV. 



CRs are reported, the first one related to subthreshold 
oscillations occurring at very low noise fluctuations in 
absence of spiking and the second one due to the reg- 
ularization of the emitted spike trains. In the regime 
of bistability the intermittent dynamics between the two 
stable dynamical states is interpreted in terms of jumping 
processes across activation barriers induced by noise fluc- 
tuations. Finally, evidences of incoherence maximization 
in the repetitive firing regime are presented. 

The model and the various employed indicators, 
namely the distribution of the interspike-interval (ISI), 
the spike-triggered averages, the autocorrelation function 
of the signal and of the ISIs, and the conditional entropies 
are introduced in Section [H] The results for the silent 
regime are reported in Section lTTTl while the bistable and 
repetitive firing regimes are examined in Section ITVl and 
fVl respectively. The last Section contains a summary of 
the results and concluding remarks. 



II. MODELS AND TOOLS 

A. The Hodgkin-Huxley Model 

The Hodgkin-Huxley (HH) model describes the dy- 
namical evolution of the membrane potential V(t) during 
a voltage clamp experiment. It can be written as 

- C^- = g Na m 3 h(V - E Na )+g K n\V - E K ) MX) 
at 

+ 9l(V - E L ) - 1(f) , 

where I(t) is an external current and the evolution of the 
gating variables X = m, n, h is ruled by three ODE's of 
the form 



X)-f3 x {V)X 



(2) 



The parameters entering in Eq. Q are C = 1/xF/cm 2 , 
E Na = 50mV, E K =-77 mV, T L =-54.4 mV, 5^=120 
mS/cm 2 , 5^=36 mS/cm 2 , and </l=0.3 mS/cm 2 . The 
expressions of the nonlinear functions ax{V) and /3x{V) 
are explicitly reported in Tabled 

We consider the single HH model subjected to Ne 
(resp. Ni) uncorrelated trains of excitatory (resp. in- 
hibitory) post-synaptic potentials (EPSPs, resp. IPSPs). 
Each post-synaptic potential (PSP) is schematized as an 
instantaneous variation of the membrane potential by a 



positive (resp. negative) amount AV for excitatory (resp. 
inhibitory) synapses. Similarly to what has been done in 
Ref. , the amplitude of each voltage kick is assumed to 
be 0.5 mV, i.e. reasonably small (« 7 %) with respect 
to the distance between the "threshold" for spike initia- 
tion for rapid EPSPs and the resting potential (~ 6 — 7 
mV) Moreover, amplitudes ~ 0.5 mV are compa- 

rable with average EPSPs experimentally measured for 
pyramidal neurons in the visual cortex of rats Q. This 
amounts to exciting the neuron with an impulsive 
current 

N e Ni 

m = o[EE^-4)-EE^-c)] (3) 

k=l I m=l n 

where t\ (resp. i^J are the arrival times of the excitatory 
(resp. inhibitory) PSPs and Q = CAV is the charge as- 
sociated to each kick. The dynamics has been integrated 
by employing a fourth order Runge-Kutta scheme with a 
time step of St — 10~ 5 — 10~ 2 ms. A spike is identified 
when V(t) overcomes a fixed detection threshold 6 = — 5 
mV. The results reported in this paper refer to averages 
performed over time spans corresponding to 30, 000 to 
600, 000 emitted spikes or (in the low noise limit, where 
spikes are more rare) at least to an integration time of 
t re 1,000- 10,000 s. 

In order to reproduce realistic inputs received by cor- 
tical neurons, for each afferent synapse the time inter- 
val distribution between PSP inputs is chosen Poisso- 
nian with an average frequency vq = 100 Hz 0. Since 
the trains coming from different neurons are assumed to 
be uncorrelated, this amounts to considering only two 
Poissonian distributed input trains of kicks of amplitude 
Ay, one for the excitatory and one for the inhibitory 
neurons with frequencies ve = N e x vq and vi = Ni x vq, 
respectively. We consider a number of input neurons of 
the order of ~ 100 — 1, 000, thus the HH neuron is stim- 
ulated with average frequencies ve(vj) ~ 10 4 — 10 5 Hz, 
consistent with a high- input regime [6j. 

The stochastic input can be characterized in terms of 
the net spike count within a temporal window AT 



N e 



7V(AT) = ]Tnf(AT) 



fc=i 



m— 1 



(AT) 



(4) 



where nf'(AT) (resp. n^(AT)) represents the number 
of afferent EPSPs (resp. IPSPs) received from neu- 
ron k (resp. to) in the interval AT. According to 
the theory of renewal processes |l| each variable n(AT) 
has a gaussian distribution with average < rt(AT) >= 
AT/aisi and variance < ?i 2 (AT) > - < n(AT) > 2 = 
(AT visi)/o,jgj, where ajsi and vjsi indicate average 
and variance of the ISI-distribution. In particular for 
Poissonian distributed ISIs with average frequency i>o' 
aisi = l/^o and v IS i = 

By assuming statistically independent input trains, 
also N(A{T)) follows a gaussian distribution with av- 
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erage and variance given by 
< N(AT) > = (N e - Ni) 



AT 
aisi 



(5) 



Var[N(AT)} = (N e + N) 



(AT v IS i) 



l ISI 



j AT 
aisi 



Here it has also been assumed that the excitatory and 
inhibitory inputs are characterized by the same aisi and 
visi- The parameter a = (y/ (N e + Nijv7si)/o>iSi mea- 
sures the standard deviation of the stochastic input pro- 
cess. 

In the high- input regime the distribution of N(A(T)) 
turns out to be the same for input ISIs following a Pois- 
sonian or a uniform distribution 



Pu(t) 



in 




otherwise 



(i+g) 



if i±ziZ < t < 

vo — — va 



(6) 



provided that both distributions have the same average 
and variance. For Pu(t), with < e < 1, ajsi = l/^Oi 
U/5j = e 2 /(fg3) and cr = Ey/(N e +iVj)/3. Since we are 
interested in analyzing the effect of the variance of the 
stochastic process on the response of the HH model for 
fixed < N(AT) >, we will consider uniformly distributed 
ISIs, because this kind of distribution allows to reach 
extremely small values of the variance (by lowering the 
e-parameter) without modifying the average. 

Both for Poissonian and uniform distributions with 
fixed aisi = l/vo the average current stimulating the 
neuron is given by 



I = 



CAV < N(AT) > 
AT 



= CAVv (N e - Ni). 



(7) 



The bifurcation parameter I determines in which dynam- 
ical regime the neuron is operating. 

Since it is equivalent to stimulate a neuron with a con- 
stant current I(t) = Idc or with a periodic train of kicks 
of sufficiently high frequency D, the frequency-current 
response curve (usually obtained with a constant input 
current) can also be recovered by considering a periodic 
input originating from Ni inhibitory and N e excitatory 
synapses each firing with a frequency vq (cf. Fig. 

Three principal dynamical regimes can be singled out 
for the HH model subject to a constant current stim- 
ulation [24|. For small currents the HH model is in 
a silent regime, i.e. its dynamics is always attracted 
by a stable fixed point. In particular, the relaxation 
towards this state is characterized by damped oscilla- 
tions, since it is a focus. By increasing the current the 
fixed point looses its stability via a sub-critical Hopf bi- 
furcation for I = Ihb — 9. 78/i A/cm 2 . The unstable 
limit cycles emerging at Ihb annihilates via a saddle- 
node bifurcation with stable periodic oscillations, corre- 
sponding to tonic firing, at a lower value of the current 
I = Isn — 6.27/iA/cm 2 . Therefore the following three 
dynamical regimes can be identified: 

1. a silent regime for / < Isn', 




I (nA/cm ) 



FIG. 1: The filled circles represent the frequency-current re- 
sponse curve (yc versus 7) obtained by considering a peri- 
odic train of kicks of amplitude AV = 0.5 mV and frequency 
v = I /(CAV) = v (N e - Ni). The solid (resp. dash-dotted) 
line indicates the stable (resp. unstable) fixed point solution. 
In the inset the time evolution of the membrane potential 
V(t), corresponding to a typical limit cycle solution, is shown. 



2. a bistable regime, where the fixed point coexists 
with a stable limit cycle solution for Isn < I < 
Ihb', 



3. a periodic firing regime for / > Ihb- 



B. Statistical and dynamical indicators 

In order to characterize the output of the neuron and 
to examine the coherence effects in the response we have 
employed the following indicators: 

• the distribution of the output interspike intervals 
Pisi(t) and its first moments: the average ISI 
(■Aisi ) and the corresponding standard deviation 

(Sisi); 

• the spike-triggered average potential (STAP) j2f| 
that gives the average shape of the membrane po- 
tential preceding the emission of a spike and the 
spike-triggered average input fluctuations (STAF): 



«(*) 



(N t (AT)— < TV (AT) >) 
y/Var[N(AT)] 



(8) 



where N t (AT) is the net spike count at time t < 
before a spike emission occurring at t = 0, while 
the expressions of < N(AT) > and Var[N(AT)} 
are reported in Eqs. JSJ. The STAF © can be re- 
lated to current fluctuations, estimated over a time 
window AT, with respect to their average value I 
(see Eq. {7J|). A positive (resp. negative) value 
of q at a certain time t < indicates a correlation 
between a positive (resp. negative) current fluctu- 
ation at that time and the emission of a spike at 
t = 0. 
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• the coefficient of variation of the ISIs 
Sisi 



R = 



A 



ISI 



(9) 



typically employed to characterize the nature of a 
process, being R = for a perfectly periodic re- 
sponse and R = 1 for Poissonian output; 

• the normalized autocorrelation function C(t) for 
the membrane potential and the correlation 
time UM defined as 



C 2 (t) dt 



(10) 



As a further indicator we have employed the condi- 
tional entropies h(N) [2^1 . In order to define these quan- 
tities, we need to digitize the output of the neuron in 
a binary sequence, where 1 and indicates respectively 
presence or absence of a spike in a certain time window 
At. The choice of the resolution At employed to analyze 
the output is crucial (2^ 1 - 

Indicating with Sk the binary symbol associated to 
the /c-th window of the time series and with Wn = 
(sj+i, Sj+2, • • • > S i+N) a sequence of symbols (word) of 
length N, the probability that a certain word will be 
observed is P(Wn)- The Shannon block entropy is then 
defined as 



H(N) = - P ( W n) l°g 2 P(W N ) 

{Wn} 



(11) 



where the sum is extended to all the 2 N possible words 
of length N. The conditional entropies are given by the 
following difference [23 



h(N) = H(N + 1) - H(N) = 

-( £ P(s N+1 \W N )\og 2 P(s N+1 \W N ) 

\{sjv+l} / 



(12) 



{W N } 



where the brackets indicate the average over all possible 
sequences of length N preceding the symbol sn+i and 
P(sn+i\Wn) is the conditional probability to observe 
symbol sjv+i once the word Wn has been registered. By 
definition we set h(0) = H(l). 



1 The time window At should be chosen in an appropriate way to 
avoid on one side to observe two spikes within the same window 
and on the other side to have extremely long series of zeros, that 
would spoil the statistical analysis. A good choice of At is the 
refractory period of the neuron following the emission of a spike 
and an upper bound for this time can be given by the measured 
minimal ISI value. Since the period of relative refractoriness de- 
creases upon increasing the amount of noise, for each current we 
have chosen At to be the minimal ISI measured for the maximum 
noise variance considered. As a matter of fact At = 5 ms has 
been used for all current values, apart from / = where At = 7 
ms was more appropriate. 



The conditional entropies represent the average infor- 
mation gained by the knowledge of the (N+ l)-th symbol 
of a time series, once the other N symbols are already 
known. For sufficiently long words, if the examined pro- 
cess can be considered as "ergodic", h(N) tends to an 
asymptotic value ft. a, the length M for which the satu- 
ration is attained gives the order (the memory) of the 
Markovian process able to reproduce the considered dy- 
namics (3 • F° r a regular process (e.g., a periodic state) 
Ha = 0, while for a purely stochastic process the condi- 
tional entropies attain their maximal value, i.e. Ha = 1 
bit 2 . 

As already mentioned, the conditional entropies have 
been successfully employed to characterize SR [2^. In 
particular a clear minimum for Ha was observed in cor- 
respondence to an optimal noise amplitude giving rise to 
a maximum in the signal-to- noise ratio (SNR), that is a 
typical signature of SR. The minimum in \ia was asso- 
ciated with the most ordered (coherent) structure of the 
output binary sequence. 

Finally, in order to characterize the correlations 
present in the binary sequence we use the following au- 
tocorrelation function 



icq 1 
2] 



_ J2i (< SjSj+k > - < Si >) 



E< (< 



< Si > 5 



(13) 



and the associated correlation time Tfo„ defined as in Eq. 

In the next sections we investigate the response of the 
HH model stimulated by stochastic spike trains in the 
three dynamical regimes shown in Fig. ^ 



III. SILENT REGIME 

In the regime where the neuron subjected to a con- 
stant current does not fire, i.e. in the range I < Isn — 
6.27/iA/cm 2 , the presence of noise in the input (charac- 
terized by the standard deviation of the noise a) induces 
stochastic firing of the model with an average firing rate 
v out = I/Ajsi steadily increasing with a (as shown in 
Fig. Hfor I = 5/^A/cm 2 ). 

For a fixed average input current / the ISI-distribution 
depends strongly on the standard deviation a of the noise 
(cf. Fig. |3J). For low noise Pisi{t) exhibits a multimodal 
structure with an exponential tail, for increasing noise 
the additional peaks and the tail disappear gradually. In 
the high noise limit Pisi{t) tends to an inverse Gaus- 
sian [jm. 



2 From a numerical point of view it is known that the estimation 
of Shannon entropy 1111 from finite samples can lead to under- 
estimation of H(N). In order to reduce the systematic errors 
arising from naive estimations of P(Wn) we have emp loye d the 
analytic estimator recently introduced by Grassberger |2SH . 
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FIG. 2: Average firing rate function of a for I = 

5/xA/cm 2 . In the inset the frequency is plotted in a lin-log 
scale as a function of 1 /a 2 for a 2 < 500. The dashed line is an 
exponential fit y = A * exp [—a 2 /a 2 ] in the range 6 < a 2 < 20 
with A = 162 and a 2 w 93. 

In the following sub-section IIII A| we discuss the mech- 
anisms responsible for the firing activity of the neuron in 
the low and high noise limits, while sub-section IIII Bl is 
devoted to coherence resonance phenomena observed for 
intermediate level of noise fluctuations. 



A. Response of the neuron to low and high noise 

Let us discuss the origin of the main features present 
in Pisi(t) at low noise. The coexistence of multi-peaks 
and the exponential tail in Fig. |2K suggests that at least 
two different mechanisms are responsible for the firing of 
the neuron. 
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FIG. 3: Probability density distribution Pisi(t) for I — 
6.15/iA/cm 2 . The fi gures refer to three different values of 
a: (a) 4.6, (b) 12.3, and (c) 100.0. The dashed line in (a) in- 
dicates the slope of the exponential tail, while in (c) it repre- 
sents an inverse Gaussian distribution with the same average 
and the same variance as the original Pisr(t). 



1. Spikes triggered by relaxation oscillations 

We can safely affirm that the multiple peaks in the 
Pjsi(t) are associated with the relaxation oscillations of 
the membrane potential towards the rest state follow- 
ing the emission of a spike. However, a linear stability 
analysis around the stable fixed point can provide the os- 
cillation periods around the focus, but is not sufficient to 
fully characterize the relaxation dynamics, that exhibits 
also nonlinear aspects. 

In order to understand if these oscillations are respon- 
sible for the multimodal structure of Pisi(t), we have 
stimulated the model (11121) with a step current of ampli- 
tude I, i.e. 



I(t) 



if t < 

1 if t > 



(14) 



and registered the shape of the output V(t). As shown 
in Fig. 0Ji (dashed line), the HH model responds emit- 
ting one or more action potentials followed by damped 
oscillations. The first oscillation of period T n i has a clear 



nonlinear character, while we have verified that the angu- 
lar frequency 1-njTi describing the subsequent oscillations 
corresponds to the imaginary part of the complex conju- 
gate Floquet eigenvalues associated to the stable fixed 
point. Therefore the latter solutions can be completely 
characterized within a linear stability analysis of model 

To better compare the periods of the relaxation os- 
cillations with the measured positions of the first (resp. 

second) peak 4sj ( res P- 4sj) °f tue ISI distribution, 
we have registered the temporal interval separating the 
peak of the action potential and the first subsequent 
maximum T n i and corrected (increased) these values by 
the corresponding rise time needed to reach the detec- 
tion threshold O. As shown in Fig. [S] the comparison 
of tfgj with the corrected values T n i is very good for 
< I < 5[iA/cm 2 . Moreover, by approaching I$n we 
observe that T n \ — * l/u c (where v c is the frequency of 
the spike train limit cycles emerging at the saddle-node 
transition). For completeness we have also verified that 
the other peaks observed in the distribution Pisi(t) are 
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FIG. 5: Positions of the peaks of the distribution Pisi(t) as 
a function of the average current I for low noise: first peak 



(i) 



(filled diamonds) and second peak tjgj (asterisks). For 
comparison the time intervals separating the induced spike 
from the maximum associated to the first (nonlinear) oscilla- 
tion T n i (empty squares) and to the second (linear) oscillation 
T n i + Ti (empty triangles) are shown. The periods \jv c of 
the regular spike trains emerging above Isn are also plotted 
(dashed line). The data for t^gj and t^j have been obtained 
for a « 2.6 - 10. 
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FIG. 4: (a) Pisi(t) obtained by the stochastic input with 
a — 4.6 (continuous curve) and the potential output V(t) 
following a step current stimulation ill Iti (dashed line). The 
position of the spike has been shifted to t = and the ac- 
tion potential has been rescaled to better reveal the relax- 
ation oscillations, (b) tjg, as a function of a. The error bars 
correspond to the histogram resolution employed to estimate 
Pisiit). Data refer to I = 6.15^A/cm 2 . 



ance, by increasing a we expect that the position of 
the first peak will be shifted to smaller values. Due to 
stronger noise fluctuations, after the emission of a spike 
a second spike can be elicited even before that the subse- 
quent relaxation oscillation of V(t) reaches its maximum 



value. This tendency of 
trated in Fig. 01). 



to decrease with a is illus- 



2. Firing activated by noise 



related to the subsequent linear oscillations of period 7] 
shown in Fig. 01i (dashed line). In particular, we con- 

(2) 

centrated on the second peak tjgj and compared it with 
T n i + Ti , also this time the agreement is quite good (see 
Fig. |SJ 3 . In general for the m-th peak we expect that 
a good approximation of its position will be given by 
T n i + (m — and this is indeed verified as shown in 
Fig. Hi. 

The damped oscillations following a spike induce a 
modulation in the degree of excitability of the model, or 
analogously a modulated effective threshold for spike elic- 
itation. Indeed a Pisi(t) quite similar to that displayed 
in Fig. QJi was obtained by Wilbur & Rinzel by consid- 
ering a leaky integrate-and-fire model with a threshold 
evolving dynamically in time p^H . 

So far we limited ourselves to low values of noise vari- 



In the first part of this subsection our analysis has been 
devoted to relatively short ISIs, however the exponential 
tail present in Pjsj (t) (see Fig. |3K) is extremely relevant 
since it gives the main contribution to the average firing 
rate v out of the neuron. As already mentioned, this long 
lasting tail should be related to some sort of activation 
process. 

An activated firing can arise in this regime due to the 
competition between the tendency of the HH dynamics 
to relax towards its stable fixed point (V res t(I)) and noise 
fluctuations that instead lead the system towards an ex- 
citation threshold (note, however, that there is no fixed 
threshold in the HH- model Q). Therefore, the dynamics 
of V(t) resembles the overdamped dynamics of a par- 
ticle in a potential well under the influence of thermal 
(stochastic) fluctuations . Due to the activation pro- 
cess the membrane potential can be driven towards the 
excitation threshold with an average escape time given 
by Kramers expression |3^ | 



3 Above the saddle- node transition I > Isn the position of the 
second peak has been estimated as l/u c + T; 



T e oc e Ws ^ 2 



(15) 



7 



where Ws plays the role of an energy barrier and a 1 of an 
effective temperature of the bath. This behavior is indeed 
verified for small a- values as shown in the inset of Fig. 
Therefore we expect that for a < a c = \/Ws (c c ~ 9.64 
for / = 5/iA/cm 2 ) the dynamics can be characterized as 
an activation process, while for a 2 larger than the barrier 
value the dynamics should be mainly diffusive. A further 
indication that for low noise the dynamics is essentially 
Poissonian is given by the fact that R « f for a < a c . 
We have estimated Ws for various values of the average 
current for I < Ihb and a linear decrease of the barrier 
height with I is clearly observable in Fig. except in 
direct proximity of Ihb- Moreover (as expected) Ws — > 
approaching Ihb where the fixed point looses its stability 
via a sub-critical Hopf bifurcation 4 . 
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emitting a spike) with higher probability when excited 
with a specific input frequency (ss 1/TJ = 61 — 88 Hz 
for < I < Isn)- This result agrees with a previous 
analysis reported in [2l| , where it has been shown that a 
silent HH neuron, subjected to a sinusoidal current, op- 
timally resonates when forced with a frequency linearly 
correlated with 1 /TJ . 
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FIG. 7: STAF q(t) (dashed line) and STAP V(t) (solid line) 
preceding the emission of a spike. At time t = the potential 
overcomes the detection threshold O. Both quantities have 
been defined by employing a time window AT = 0.1 ms and 
have been averaged only over ISIs longer than 150 ms. The 
data refer to I = 5fiA/cm 2 and a — 5.7. For a better com- 
parison q(t) has been shifted and amplified by a factor 10, the 
original function is displayed in the inset. 



FIG. 6: Activation barrier heights W as a function of the 
average input current for the silent state Ws (filled circles) 
and the oscillatory state Wo (empty squares). The vertical 
dashed (resp. dot-dashed) line indicates the position of 7s jv 
(resp. Ihb)- A linear fitting to Ws in the interval 2fi A/cm 2 < 
I < 9 //A/cm 2 is also shown as a dotted line. 

To gain some deeper insight into the role of noise fluc- 
tuations in eliciting a spike, we have estimated the STAP 
V(t) and the STAF q(t) for sufficiently long ISIs and for 
small a. As shown in Fig. [7| q(t) and V(t) exhibit os- 
cillations of period « T; preceding the spike emission 
and these oscillations are almost in phase, apart from 
a small delay in the response of V(t) with respect to 
the input noise. This indicates that the potential fol- 
lows the current oscillations and that the emission of a 
spike (for long ISIs) is commonly triggered by the exci- 
tation of linear subthreshold oscillations around the rest 
potential. Therefore the neuron in proximity of the rest 
state acts as a sort of selective filter since it responds (by 



4 The average escape times T e have been estimated directly as 
for I < Isn since at sufficiently low noise the activation 
is the prevailing mechanism in this regime; above Isn where 
there is a coexistence of a stable limit cycle with a stable fixed 
point these times have been measured as the residence times Ts 
in the silent state (i.e. in proximity of the fixed point). 



3. High noise limit 

As already mentioned, for sufficiently high noise the 
probability density distribution of the ISIs reduces to an 
inverse Gaussian. The reason is that for a a c the 
"potential well" surrounding the stable fixed point be- 
comes irrelevant for the dynamics due to the amplitude 
of the noise fluctuations, whose effect on the dynamics 
of the neuron are twofold: a constant current I driving 
the system plus a stochastic (Gaussian) term with zero 
average. Therefore we can assume that the dynamics 
of the membrane potential is determined by a Langevin 
equation with a drift term. The distribution of the first 
passage times from a certain threshold, in presence of a 
reset mechanism to a fixed rest potential, is give n by the 
well known inverse Gaussian distribution [lLl30| 



/(*) 



2/3t 



(16) 



where a — Ajsi and (3 — Sj SI /Ajsi- The good agree- 
ment shown in Fig. (5J; confirms that the mechanism 
leading to repetitive firing in the high noise limit can 
be indeed schematized as a Wiener process plus drift. 

Moreover, by assuming that the drift is proportional 
to / and the amplitude of the effective noise to a one 
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obtains the following dependence for the ISI coefficient 
of variation 



R 



(17) 



(7 + 7 )V^7 ' 

where Iq is a parameter. For fixed noise variance, R x 
y/Aisi should then be inversely proportional to I + Iq as 
indeed verified for a — 100 in Fig. EI 




I (uA/cm ) 



FIG. 8: R x v^47s7 (filled circles) versus I, the dashed line 
is a best fit to the data of the form a/(I + lo) with a = 
16.35 ms 1 ' 2 [lA/cm 2 and Jo = 14.59 /jA/cm 2 . The data refer 
to cr = 100. 



B. Coexisting coherence resonances 

In the silent regime we have encountered two coexisting 
coherence resonance phenomena. The first one is related 
to the existence of an optimal noise level for the regu- 
larization of the output spike trains and it corresponds 
to the effect previously reported in pol l2lj . The sec- 
ond coherence effect can only be detected by considering 
the membrane potential dynamics, not being associated 
with spikes, but instead with the excitation of quite reg- 
ular sequences of subthreshold oscillations occurring at 
small noise variances. 



1. Coherence of the emitted spike trains 

At intermediate noise levels the Pisi{t) reduces its ex- 
ponential tail and begins to assume the shape of an in- 
verse Gaussian distribution (see Fig. |2b). Therefore, the 
activation mechanism responsible for the firing is gradu- 
ally substituted by another stochastic mechanism resem- 
bling a Wiener process with drift. In correspondence to 
the transition from one kind of stochastic process to the 
other a regularization of the output signal is observed: 
this phenomenon is known as coherence resonance [l3| 
and has already been reported for various neuronal mod- 
els. Our aim in the present subsection is to give a more 



detailed characterization of the CR phenomenon in terms 
of commonly used indicators of coherence [l4|, like the 
coefficient of variation R and the correlation time r c , but 
also in terms of the saturated conditional entropies Ha, 
not previously employed in this context. 



a=5.7 
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FIG. 9: Membrane potentials V(t) for various levels of noise 
a in the silent regime at I — b^A/cm 2 . 

Let us first examine the signals displayed in Fig. [5] 
for I = 5fiA/cm 2 . At low cr-values one observes rare 
spikes induced by the activation mechanism, while for 
increasing noise the train becomes more and more regular 
and finally at high cr-values the noise begins to modify 
even the relative refractory times. 

A first quantitative evidence of coherence resonance 
can be given in terms of R: at low a this quantity tends 
to one as expected for Poissonian processes, indeed R 
reaches a value slightly greater than one due to the mul- 
timodal structure of the Pisi(t); for increasing a the SVsj 
decreases faster than Aisi since the diffusive nature of 
the process begins to prevail over the activation effects 
(exemplified by a Poissonian distribution for Pisi(t) with 
Si si — Aisi). The behavior of R at large cr-values is 
given by the expression (|17fl and since Aisi slightly de- 
creases for large noise amplitude the coefficient of vari- 
ation turns out to be an increasing function of a. A 
minimum of R is found at intermediate noise levels due 
to the predominance of two different mechanisms at the 
origin of spiking in two opposite limits: activated barrier 
crossing in the limit a << cr c and diffusive motion with 
drift in the opposite limit a >> cr c . Indeed, as shown 
in Fig. 1101 a minimum is attained for a = <tr ~ 55 for 
I = 5[iA/cm 2 . 

A second striking evidence of coherence resonance is 
given by a maximum in the correlation time r c estimated 
directly by integration of the squared autocorrelation 
function of the signal according to Eq. (|10fl . This max- 
imum occurs for a value smaller than that obtained for 
R, i.e. a = a T w 30 for I = 5fiA/cm 2 . 

As previously suggested for SR [2^], dynamical en- 
tropies are indicators appropriate to characterize the de- 
gree of unpredictability of the output of a certain non- 
linear system, therefore we expect that they should be 
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FIG. 10: / = 5/xA/cm 2 : coefficient of variation R (empty 
circles) and correlation time r c (filled circles) as a function of 



useful also in the present context. In Fig. ^] the con- 
ditional entropies h(N) are reported for increasing word 
length N at I = 5/iA/cm 2 . These quantities saturate to 
their asymptotical values Ka for N — 5, thus indicating 
that the present dynamics can be reproduced in terms of 
a Markovian process of order 5 with a temporal memory 
Rj 5x At = 25ms. Since for this current value Aisi varies 
between 17 ms and 11 ms for 20 < a < 150, this means 
that in order to recover all the transmitted information it 
is sufficient to record from two to three successive spikes. 
Similar to what has been reported in Ha exhibits a 
non-monotonic behavior characterized by a rapid initial 
growth at small a followed by a decrease and a minimum 
at intermediate noise (a = ah ~ 33). Initially the signal 
is extremely regular, being essentially a sequence of zeros 
in absence of spikes, then the noise tends to randomize 
the system and this induces a rapid increase. Finally, 
the signal regularizes at finite a and this leads to the 
occurrence of a minimum in the asymptotic entropies. 




FIG. 11: Conditional entropies h(N) as a function of a for 
various word lengths TV. The data have been obtained for 
I — 5fj,A/cm 2 with a resolution of At = 5 ms. 

We can safely affirm tha the three employed indicators 



agree in indicating a clear coherence resonance effect at 
the considered current / = 5fiA/cm 2 . Furthermore, we 
performed a similar analysis at various current / in the 
silent and in the bistable regime and the corresponding 
values <tr, ah, a T are shown in Fig. 1121 For all indica- 
tors the optimal amount of noise needed to observe CR 
decreases for increasing current. However, while ctr and 
ah seem not to vary dramatically in the examined range, 
a T decreases noticeably approaching the SN bifurcation. 
This behavior seems to be due to the fact that the signal 
autocorrelation function does not only register coherence 
effects associated to the ISIs, but also those related to 
the subthreshold oscillations occurring at quite low noise 
(as reported in the next subsection). The interaction 
between these two resonances leads to the enhanced de- 
crease of a T , somehow the ISI resonance is entrained by 
the second one and shifted towards small a. Other indi- 
cations that the signal becomes more and more correlated 
approaching Is n are given by the corresponding increase 
of the maximal correlation time and by the decrease of 
the minimum of R and of Iia (all measured at the res- 
onance). This scenario is consistent with the fact that 
the activation barrier Ws separating the rest state from 
the excited state also decreases with I and therefore the 
latter becomes more accessible at lower noise variances. 




I (u A/cm 



FIG. 12: Optimal noise standard deviations corresponding to 
coherence resonance versus I for the three considered indica- 
tors: namely, (filled circles), o> (filled triangles) and o~h 
(empty squares). The vertical dotted line indicates Isn, while 
the dashed lines are guides for the eyes. 



2. Coherence of the subthreshold oscillations 

Let us now discuss in more detail the origin of the 
second coherence resonance observed at very low noise 
variance. In Fig. ^| the behavior of the correlation time 
t c as a function of a is shown for / = 4/iA/cm 2 in a wider 
range of noise with respect to the data shown in Fig. ITU1 
In that figure the examined noise range was restricted 
to values for which the statistics of emitted spikes was 
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sufficiently rich to ensure a meaningful definition of R. 
In Fig. 113b r c reveals two clear maxima: the higher one 
located at a « 3 and the lower one at a ~ 33. The origin 
of the first peak can be understood as follows: for a < 3 
almost no spikes are emitted by the neuron, however, the 
increase of noise tends to stimulate series of subthresh- 
old oscillations that are more and more correlated; for 
a > 3 the statistics of the emitted spikes begins to be no 
more negligible and the occurrence of rare spikes tends 
to decorrelate the signal leading to a decrease of r c . By 
further increasing the noise variance the signal begins 
to be characterized by sequences of spikes, therefore the 
autocorrelation function starts to register essentially the 
correlation of these events and r c reveals a second peak 
related to the regularization of the spike trains, this is 
the coherence phenomenon previously discussed. 

The transition from one dynamical regime to the other 
can be better understood by examining the autocorrela- 
tion function of V(t) for various values of a. As shown 
in Fig. 113b . for a ss 3 the autocorrelation function C(t) 
reveals oscillations of period w TJ, while at a w 9.7 the 
maxima of C (t) are located in correspondence with mul- 
tiples of t « tj)j. In between these two values there is a 
transition from an output signal dominated by the sub- 
threshold oscillations to a signal characterized by trains 
of spikes, this transition can be located at a w 4.5, since 
at this value C(t) exhibits, at the first oscillation, two 
maxima one located at t ~ Tj and one at t sa tjgj. 

To summarize, for the HH model in the silent regime 
two kind of coherence resonances can be observed, one 
at quite low noise related to the excitation of subthresh- 
old oscillations around the rest state and another one at 
higher noise due to the regularization of successive I Sis 
associated to the spikes emitted by the neuron. Obvi- 
ously, since the first resonance is not related to spike 
occurrence it cannot be revealed by R or by h, but only 
by t c . 



IV. BISTABLE REGIME 

In the bistable regime the modifications of the shape 
of the ISI distributions due to the effect of noise resemble 
those found in the silent regime, apart for very low noise 
variances (see Fig. 1140 . As shown in Fig. 114b . in this 
latter case one observes a quite pronounced peak corre- 
sponding to l/v c followed by smaller multi-peaks and an 
exponential tail. This reflects the fact that the dynam- 
ics due to the noise switches back and forth between the 
two coexistent states: oscillatory and silent, (cf. the 
time evolution of the membrane potential reported in 
Fig. ED in particular at a = 3.2 (for I — 8/xA/cm 2 )). 
This dynamical behavior is reminiscent of the motion of 
one particle in a double well subjected to thermal fluctu- 
ations. Therefore, the residence times in the two states 
should be related with the noise amplitude by a Kramers 
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FIG. 13: I = 4fiA/cm 2 : (a) r c (filled circles) as a function 
of a, the solid line is a guide for the eyes; (b) autocorrelation 
functions C(t) for three different noise values (namely a = 
2.9,4.5 and 9.7), the vertical dotted line indicates the period 
Ti, while the vertical dashed line refers to for a — 9.7. 



relationship in the standard manner 

T.oce^ 2 , (18) 

where x = S for the silent state and x — O for the os- 
cillatory state. Moreover, the distributions of Ts (resp. 
To) should be Poissonian due to the stochastic nature 
of the jumping from one state to the other. In order 
to verify Eq. 1)18(1 we have estimated the times Ts and 
To as a function of a for small values of this param- 
eter. The data, displayed in Fig. ^3 confirm that in 
the low noise limit both times can be described in terms 
of an activation process induced by noise fluctuations. 
Moreover, barriers Wo and Ws have been estimated 
for various currents and are reported in Fig. For 
Isn < I < 8/zA/cm 2 the most stable state is the fixed 
point, while above / = 8fiA/cm 2 the stability of the os- 
cillatory state prevails. 

Let us now come back to the Pisi(t) for a = 2.2, in this 
case we have estimated the distributions of the ISIs dur- 
ing the oscillatory state only, let us term this probability 
distribution as Pjgj{t). This distribution, reported in 
Fig. 114b . reveals a multimodal structure. While the first 
peak corresponds to « X/v c , the other peaks are related 
to the linear damped oscillations towards the fixed point, 
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FIG. 14: Pisiit) in the bistable regime at / = SfiA/cm 2 . 
Pisi(t) is displayed for various values of the noise: a = 2.2 
(a), a = 7.3 (c), and a = 30 (d). Additionally in (b) P r ISI {t) 
is shown for a = 2.2. 




FIG. 16: Residence times To and Ts as a function of 1/cr 2 
for low noise variance values: the dashed lines are fits to the 
numerical data with expression 11811 for a < 4, the corre- 
sponding barriers are Ws = 21.4 ± 0.5 and Wo = 17.9 ± 0.3. 
All data refer to / = 8fiA/cm 2 . 
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FIG. 15: Membrane potential V(t) evolution in time for var- 
ious value of the noise variance. The top signal refers to 
I = 8fj,A/cm 2 and a = 3.2, the other three curves corre- 
spond (from top to bottom) to a = 6.8, 38.7 and 141.6 and 
/ = 7/i A /cm 2 . 



i.e. tjt^j 1 ^ — tfgj ~ kTi for k > 1. The origin of this mul- 
timodal structure can be understood by comparing the 
STAF qi (t) preceding the emission of a spike for interme- 
diate ISI-values (namely for ISI durations corresponding 
to the second and third peak of the distribution) with the 
STAF qh{t) associated with long ISIs falling in the expo- 
nential time tail. Analogously to what observed in the 
silent regime qh{t) reveals oscillations of period « T;, in- 
dicating that the neuron will return to fire, once entered 
in the silent state, mainly when stimulated via current 
oscillations of the proper period and in phase with V(t) 
oscillations (see Fig. 0). As shown in Fig. El Qi(t) has 



a quite peculiar behavior, it also exhibits oscillations of 
period i=s 7} , but these oscillations are in anti-phase with 
respect to those of <?l(£), apart from the last two pre- 
ceding the firing of the neuron. This means that, once 
entered in the oscillatory state, the neuron can become 
silent only if a series of current oscillations (in anti-phase 
with respect to the damped oscillations of V(t) follow- 
ing a spike) inhibit repetitively spike emission allowing a 
sufficient decrease of V(t) leading the system in the at- 
traction basin of the rest state. If by chance one of these 
oscillations returns in phase with V(t) (during the relax- 
ation) the neuron will fire immediately after (as shown in 
Fig. I17fl . This is the mechanism at the origin of the var- 
ious peaks present in Pisi{t) at low noise values, some 
indications consistent with this scenario have been re- 
cently reported in Ref. ||- 

As far as CR is concerned, a clear resonance is observ- 
able by inspecting the dependence of R and the saturated 
\ia as a function of noise intensity for / = 7fiA/cm 2 (see 
Fig. 1181) . On the other hand, t c exhibits two almost coin- 
ciding maxima at very low noise: one due to subthreshold 
oscillations at a = 1.8 and the other one related to ISI 
coherence at a T = 2.6. At higher currents the two max- 
ima merge and they are no more distinguishable one from 
the other. In order to understand if a ISI coherence effect 
is still clearly discernible at the level of correlation func- 
tions we have estimated the time Tbi n defined in terms 
of Can (|13fl . This quantity, which is not influenced by 
the subthreshold oscillations, shows a clear maximum at 
a = 20 (see Fig. IT%K). Since Chin and \ia measure both 
dynamical features of the ISI sequence we expect that 
they should give similar information, indeed for this cur- 
rent <Jh = 26 ± 2. These results seem to indicate that the 
information reduction needed to pass from the signal to 
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FIG. 17: STAF for long ISIs q L (t) (dashed line) and inter- 
mediate ISIs qi{t) (solid line) for 7 = 8/xA/cm 2 and a = 3.2. 
qh(t) has been estimated by averaging over I SI > 70 ms, 
while qi(t) refer to 25ms < 757 < 70 ms. 



the ISIs allows to better single out the "standard" coher- 
ence phenomenon in this regime. By further increasing 
the current the only indicator that continues to signal a 
coherence resonance effect is R, while at I = 9/x Aj cm 2 
the minimum present in h becomes a shoulder and also 
Tu n does not exhibit maxima at intermediate noise. The 
main difference between the dynamics at I — 7fj,A/cm 2 
and I = 9 /i A /cm 2 is related to the fact that the oscilla- 
tory state becomes more stable than the silent state at 
this latter current (i.e. Wo > Ws)- To summarize, the 
coexistence of two coherence resonance is observed also 
in the bistable regime until the fixed point remains the 
most stable solution (i.e. for I < 8fj,A/cm 2 ). 
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FIG. 18: Regime of bistability (7 = 7fiA/cm 2 ): R (a) , r c (b), 
Ka = h(5) (c) and Tbm (d) as a function of a. 



V. REPETITIVE FIRING REGIME 

Above Ihb the system is in a stable regime of repet- 
itive firing if subjected to a constant current. As we 
can see from Fig. the noise has different effect this 
time. In particular, for sufficiently small a we observe 
that Pi siif) is essentially a Gaussian centered around 
the repetitive firing period \/v c (see Fig. I19fa). Initially 
for 1 < er < 3, Si si increases linearly with cr, as expected 
for additive noise of sufficiently small intensity. Upon 
further increasing the variance of the noise, the distribu- 
tion becomes again multimodal (as shown in Fig. 119b ) 
with peaks located at integer multiples of 1/V C , but this 
time the activation tail is almost absent. The multi-peak 
structure (mainly limited to two peaks only) is due to the 
fact that sometimes the neuron fails to emit a spike at 
1/V C because of a current fluctuation in anti- phase with 
the suprathreshold oscillations of the membrane poten- 
tial. This mechanism is analogous to the one already 
described in the previous Section to explain the multi- 
peaks observed in the bistable regime. At larger a the 
second peak reduces to a shoulder (see Fig. 1191: ) and then 
Pisi(t) converges towards an inverse Gaussian. 
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FIG. 19: Pisi(t) corresponding to I = U/iA/cm 2 for various 
cr- values: a = 2.2 (a), 8.6 (b), 12.3 (c) and 22.4 (d). 

In this regime the only indicator giving evidence of 
CR is the coefficient of variation 7?, while for a — > the 
correlation time r c diverges to infinite, since the system 
converge to a stable limit cycle, and the conditional en- 
tropies do not exhibit any relative minimum. As shown 
in Fig. 1201 7? has a maximum at a = 12.3 (for 7? w 0.19) 
followed by a minimum at a = 30 (corresponding to 
7? w 0.17). However, due to the very limited variation 
of 7?, it is difficult to appreciate from the signal evolu- 
tion some difference between a = 12.3 and 30. A similar 
effect of maximal spike train incoherence was observed 
for a leaky integrate-and-fire model with an absolute re- 
fractory period for supra-threshold base current [T^ . In 
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that case a maximum occurs at intermediate noise values 
since perfectly regular spiking was found for vanishing 
noise and in the large noise limit. In our model the max- 
imum in R is related to the emergence of the multi-peak 
structure of Pisi{t). In the limit a — > we observe regu- 
lar spiking (i.e. R — ► 0), the introduction of noise in the 
system initially leads to small irregularities in the spiking 
(reflected in the linear increase of Si si with a) and then 
to the emergence of the second peak in Pisi{t), that in- 
duces an abrupt growth in the 57s/-value (as shown in 
the inset of Fig. I20|) and also in R. At higher noise levels 
(cr > 23) the dynamics reduces essentially to a Wiener 
process plus drift, this implies a merging of the two peaks 
associated with a simultaneous decrease of Sisi followed 
by a saturation. The further randomization of the dy- 
namics leads to a decrease of Aisi, that tends to its 
asymptotic value, which is the refractory time. In addi- 
tion, this behavior of Aisi is responsible for the minimum 
of R at a = 30 and for the successive growth consistent 
with Eq. l(T7j) . 
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FIG. 20: R (filled circles) versus cr for J = llfiA/cm 2 . In 
the inset are reported Aisi (empty squares) and Sisi (filled 
triangles) for the same current. 



VI. SUMMARY AND FINAL REMARKS 

In the present paper we have studied the Hodgkin- 
Huxley model in the high-input regime subjected to 
stochastic trains of uncorrelated inhibitory and excita- 
tory kicks. Our analysis suggests that the specific shape 
of the time distribution of the arrival times play no role in 
this framework. Moreover, the response of the model is 
completely determined once the average and the variance 
of the stochastic input are given. 

In the silent regime we have reported the coexistence 
of two coherence resonances: one corresponding to the 
rcgularization of the emitted spike trains at intermediate 
noise levels and the other one to the stimulation of sub- 
threshold oscillations occurring at very low noise. All the 
employed indicators, namely the coefficient of variation, 



the signal and ISI correlation times, and the conditional 
entropies, are able to identify the first resonance, while 
the second one can be detected only by the autocorrela- 
tion time of the signal. Conditional entropies, used for 
the first time in the context of coherence resonance, turns 
out to be an effective indicator for coherence. Experi- 
mental evidences of a similar coexistence of two kinds of 
resonance has recently been reported in |33j for measure- 
ments performed on an electrochemical cell and numeri- 
cally corroborated by simulations of a FitzHugh-Nagumo 
model. However, there are two differences with our anal- 
ysis: in Ref. |33| the subthreshold signal was injected 
in the system together with the noise (therefore it is a 
stochastic resonance, and not a coherence resonance) and 
the two resonances are both observable at the level of 
spike trains. 

The dynamics in the bistable regime can be described 
by activated jumping processes across barriers between 
two stable solutions : namely, the oscillatory and silent 
state. The relative stability of the two dynamical regimes 
is ruled by the ratio of the corresponding barrier heights: 
until the fixed point solution remains the most stable 
state the observed dynamics resembles that in the silent 
regime. In the repetitive firing regime the only notice- 
able feature is related to a maximization of incoherence 
observed at a finite noise level. 

Moreover, we have clarified the various mechanisms 
responsible for spike triggering and for ending repetitive 
firing. In particular, at relatively low noise the silent neu- 
ron can fire due to stochastic fluctuations via two mech- 
anisms: one related to relaxation oscillations following 
a spike and another one associated with noise induced 
activation processes. The first mechanism lead to the 
multi-peaked structure of the ISI distributions while the 
latter is responsible for the exponential tail. The pres- 
ence of peaks in the Pisiif) suggests that the system, un- 
der the influence of stochastic inputs, can resonate when 
forced with specific frequencies corresponding to 
the main peak (due to nonlinear effects) being associated 
to frequencies in the 7-range 0| (namely, from 40 to 66 
Hz for I € [0 : 9] /j,A/cm 2 ), while the second one to lower 
frequencies (namely, from 30 to 37 Hz for the same in- 
terval of currents). Indeed these results can represent an 
explanation of recent findings |35j , where clear stochastic 
resonance effects were observed in the range from 30 to 
65 Hz for a Hodgkin-Huxley neuron subjected to Pois- 
son distributed trains of EPSPs and IPSPs with peri- 
odically modulated rates plus a subthreshold harmonic 
signal. The authors [35j suggest that the origin of such 
SR should be related to subthreshold oscillations, how- 
ever the corresponding frequencies (i.e. 61 < 1/TJ < 92 
Hz for I 6 [0 : 9] /lA/cm 2 ) are too high to match with 
the observed resonance. Once relaxed in the rest state 
the Hodgkin-Huxley neuron begins to act as a selective 
filter responding to current fluctuations with the same 
frequency of the linear subthreshold oscillations, in agree- 
ment with the analysis in plj . 

In the bistable and repetitive firing regime the multi- 
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modal structure of Pisi(t) at low noise is related to the 
periodicity of suprathrcshold and (linear) subthreshold 
oscillations. Tonic firing states end when the neuron is 
stimulated with fluctuations in anti-phase with respect 
to the internal oscillations. In the high noise limit the 
neuron dynamics in all the three examined regimes can 
be represented as a stochastic Wiener process plus drift. 

In conclusion, we have shown that the Hodgkin-Huxley 
neuron in the high input regime displays a large variety of 
dynamical behaviors, thus rendering the study of its dy- 
namics interesting per se and not only for its biophysical 
implication. The richness of the Hodgkin-Huxley dynam- 
ics is particularly pronounced in the silent regime for low 
input noise variance, where the response of the model 
resembles more the activity of cortical neurons, since it 
is characterized by an almost Poissonian distribution of 
interspike intervals @|. In this regime the single neuron 
response can selectively enhance stochastic stimulations 
corresponding to different spectral ranges, thus allowing 
for a contemporary detection and transduction of a va- 
riety of signals. Moreover, a network of such neuronal 
elements will have the capability to exhibit coherent and 



correlated activity over different time scales (mainly in 
the 7 and /^-ranges 0), a property that is believed to be 
important for information encoding for cortical process- 
ing j3(|. Indeed, it has been found [2l| that a globally 
coupled Hodgkin-Huxley network subjected to stochastic 
inputs reveal, for sufficiently strong synaptic coupling, 
a collective synchronized rhythmic firing in a range of 
40 — 60 Hz, induced via coherence resonance. 

As future developments of the present work, we plan 
to investigate the role played by correlations among the 
synaptic inputs in enhancing or depressing the coherence 
effects here discussed. 
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